Universal equilibrium distribution after a small quantum quench 
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A sudden change of the Hamiltonian parameter drives a quantum system out of equilibrium. For 
a finite-size system, expectations of observables start fluctuating in time without converging to a 
precise limit. A new equilibrium state emerges only in probabilistic sense, when the probability 
distribution for the observables expectations over long times concentrate around their mean value. 
In this paper we study the full statistic of generic observables after a small quench. When the quench 
is performed around a regular (i.e. non-critical) point of the phase diagram, generic observables are 
expected to be characterized by Gaussian distribution functions ("good equilibration"). Instead, 
when quenching around a critical point a new universal double-peaked distribution function emerges 
for relevant perturbations. Our analytic predictions are numerically checked for a non-integrable 
extension of the quantum Ising model. 
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Introduction Imagine to prepare a closed quantum 
system in a given initial state po and let it evolve freely. 
After waiting a sufficiently long time an equilibrium, av- 
erage state p emerges. Because of the unitary nature 
of the dynamics, in a finite system, the evolved state 
p (t) cannot converge to p either in the strong nor in the 
weak topology [26]. Equilibration in isolated quantum 
systems only emerges in a probabilistic fashion. We say 
that the observable O equilibrates to O if the expectation 
value (O (t)) spends most of the times close to its average 
O. In other words, (O (i)) is seen as a random variable 
equipped with the (uniform) measure dt/T in the inter- 
val t £ [0, T] where T is the total observation time which 
will be sent to infinity. The probability distribution of O 
is P (o) := 5 (o — (O (t))), where the bar refers to tem- 
poral averages: / := limT^cc .T^ 1 J Q T / (t) dt. Broadly 
speaking concentration phenomena for P (o) correspond 
to quantum equilibration. The average value of a generic 
observable is readily obtained as O := (O (t)) = tr (~p~0), 
an e quat ion that defines the equilibrium state to be 
~p = p(t). Equilibration however, is related to the con- 
centration of the distribution P (o) , a convenient def- 
inition of which is encoded in the variance AO . In 
Ref. [J it has been shown that the variance of any 
observable is bounded by the purity of the equilibrium 
state V (Jo) := tr(p 2 ): This is an encouraging result, if 
V (p) is small one has equilibration for every observable. 
Equilibration should depend on the dynamic and possi- 
bly on the initial state, not on the specific observable. 

A convenient setting to probe quantum equilibration is 
that of a sudden quench. The system is initialized in the 
ground state of some Hamiltonian H\ , and then evolved 
unitarily with a small perturbation H2 — H\ +SXV. This 
situation is compelling both from a theoretical and an 
experimental point of view thanks to the recent advances 
in cold atoms technology [3, 0, HJ . 

In this paper we will analyze the full statistic of a 
generic observable P (o) after a small quench. For small 
quenches performed around a regular (i.e. non-critical 



point) the expected distribution P (o) is Gaussian in the 
generic case. Equilibration is achieved in a standard fash- 
ion. Instead for quenches performed around a critical 
point the distribution of generic observables tend to a 
new, universal double peaked function which we are able 
to compute. 

This behavior has been first demonstrated in [|| for a 
particular observable (the Loschmidt echo) on the hand 
of an exactly solvable model (Ising model in transverse 
field). Here we show that the scenario first advocated 
in [6j is in fact general to small quenches for sufficiently 
relevant perturbations. 

Critical scaling of the time-averaged state Here we 
consider the equilibrium distribution for small quench. 
When the quench is small one can either expand the 
eigenvectors of the evolution Hamiltonian H2 with per- 
turbation +6 XV or expand the initial state with respect 
to a perturbation — 8XV . We take the latter point of view. 
Let the t > Hamiltonian be H 2 = ^2 n E n \n)(n\. The 
initial state \tp ) is the ground state of Hi = H 2 — 5XV . 
Then 



hM = io)+5A]r 



n#0 



- R 



(2) ' 



O (5X 2 ) 



(note the plus sign in V). If the spectrum is non- 
degenerate the equilibrium state has the form p = 
J2 n Pn\ n )( n \ [JUL El- The weights, up to second order in 
the quench potential, are given by 
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Note that up to the same order, the purity of the equi- 
librium state is given by tr (p 2 ) = p 2 ,. The weight p is 
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precisely the square of the well studied ground state fi- 
delity F = |(0|i/) n )| @, II I! E3] and its scaling properties 
are well known [11[ . If the perturbing potential is exten- 
sive and the quench is done around a regular (i.e. non- 
critical) point F ~ exp (—const x 5X 2 L d ) where d is the 
spatial system dimension. Instead for quenches at a crit- 
ical point F ~ exp (-const x 5X 2 L 2< - d+ ^- Av ^) , ( is the 
dynamical critical exponent and Ay is the scaling dimen- 
sion of the perturbation V. Indeed it is intuitively clear 
that by shrinking 5X at will one should be able to trans- 
fer most of the spectral weight to p , a limit in which the 
purity is large. The above scalings tell us that we must 
have SX < L-Q with Q = d/2 (Q = d + ( - Ay = l/v) 
in the regular (critical) case. These are the regimes of 
small quench characterized by a large purity and hence 
large variances for generic observables. In other words 
poor equilibration. 

However the distribution of the p n for critical and reg- 
ular quenches are radically different. As we will see, this 
has direct consequences to the general form of the distri- 
bution of generic observables. 

In case of a critical quench there exist modes with 
vanishing energy: Ek — Eq = vk 1 * where k now is a 
quasi-momentum label. According to Eq. ([1]) the cor- 
responding weight pk becomes large and might even (ap- 
parently) diverge when k — + 0. In a finite system with 
periodic boundary conditions the momenta are quantized 
as fc = 27rn/L, then one would infer that, for a certain 
weight pi ~ 5X 2 L 2< *. This, however is not the correct 
scaling as we did not include the scaling of the matrix 
element. To find the exact scaling we can reason as fol- 
lows. Define the functions M (E n ) := 5X 2 \(0\V\n)\ 2 , and 
p{E n ) :— p n . With the help of the density of states 
p (E) = trS (E — H), one can write the fidelity suscepti- 
bility x as 

X = 5X 2 } p = / K —^p(E)dE. 

^ (E m -E Q f J El (E-Eof 

(2) 

We are interested in the scaling properties of M (E) after 
a rescaling of the energy. At criticality it is natural to 
assume that M (E) be an homogeneous function at the 
lower edge: M (E) ~ (E — Eo) a . Instead the product 
p (E) dE is invariant under rescaling of the en ergy . The 
scaling of the fidelity susceptibility is known \ ~ 

L 2(d+C-A v ) „ £-2(d+c-Av)/c S0; from x „ E a ~ 2 , we 

obtain a = 2 (Ay — d) /£. Using the fact that, for the 
operator driving the transition Ay = d + Q—l/v (l2j |, we 
obtain 

p{E) ~<5A 2 £T 2 /(C"), (3) 

In this last equation the energy is measured from the 
ground state, so that, being the system critical, E can 
be arbitrarily close to zero in the large size limit. The 
prediction Eq. ([3]) agrees with an explicit calculation on 
the quantum Ising model (p (uj) — 2c(ui) in Q) 

As a by-product of this analysis we obtain (0|V|fc) ~ 
6XL d - A ' ; = 8XL~ C ' +1 / U . Note that here V is the exten- 



sive perturbation. If V = ^2 X V (x), for the intensive 
component we get 

(0|V (a:) \k) ~ 5XL~ Av = 5XL-^ d+1 ' 1 ' . (4) 

Equation Q is in agreement with the analysis of [l3j per- 
formed on the sine-Gordon model. In that case d = (, = 1 
and one gets (01 cos (0<j> (x)) \k) ~ L~ 2+1 I 1 '. In fact for- 
mula (12) of [l3| can be written as (0| cos (/3(f) (x)) \k) ~ 
L~ K where K = 2- l/v is the scaling dimension of the 
cosine term. 

The content of equation © is the following. For a 
relevant perturbation (d + £ > Ay) of a critical point 
some spectral weights p n tend to be large. At finite size, 
the lowest modes have energy, E n — v (2ttti/L)^ so that 
p n ~ 8X 2 L 2 I V . In practice, since in the region of validity 
of perturbation theory, po is already "large" , the sum rule 
TlinPn = 1 constrains to have only very few p n apprecia- 
bly different from zero. We expect this scenario to be 
more pronounced for strongly relevant perturbations, in 
other words when the exponent 2/v is large. When this is 
the case, the sum rule can be saturated by taking a very 
small number of terms n max : 1 = Y, n P n ~ En=o _1 Pn- 
In our numerical simulations (see below) we have verified 
that for a case with v = 1 the sum rule is already sat- 
urated by taking as little as three terms i.e. rt max = 3. 
Moreover most of the weight is splitted between po & n d 
pi, while pi is already orders of magnitude smaller. 

The same considerations can clearly be drawn for 
the amplitudes c„ = (n|-0o) = 5A(n|V|0)/ {E — E n ) + 
O (<5A 2 ) for n > 0, for which p n = |c„| 2 . Defining the 
function c(E n ) = c n with the same reasoning as above, 
one sees that, for £J — * 0, c(E) ~ SXE- 1 /^ . Al- 
ternatively, for some low lying excitations with quasi- 
momentum k, Ck — (k\tpo) ~ SXL 1 /". Since c(E) is a 
rapidly decreasing function, and because of the sum rule 
for the Cjj, one obtains a good approximation for the time 
evolved wave- function by just resorting to very few, n max , 
amplitudes: \ip (t)) « X)n=o x_1 c n e~ %tEn \n). 

Equilibrium distribution for small quenches Let us 
now illustrate what are the consequences of these find- 
ings on the equilibration. Consider the time evolution of 
a generic observable (O (t)} . We will also give results for 
the Loschmidt echo (LE) as it is attractingan increas- 
ing amount of attention [3 EE GJ, EE El M HI • The 
Loschmidt echo is defined as C(t) — | (ijjo\e~ ltH2 \tpo) \ ■ 
Note that, as pointed out in Q the LE can be written as 
the expectation value of a particular observable (Oc (t)) 
with Oc given by Oc = IV'oKV'oI- Expanding C(t) and 
(O (i)) in the eigenbasis of H2 we obtain: 

£{t)=£+J2 2 PnPm COS [t {E n - Em)] (5) 

m>n 

0(t)=0+J2 (n\0\m)c m c^e- u ^- E ^ 

= 0+Y, 2(n|0|m)c ro c n cos [t (E m - En)] . (6) 
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Where in the last line we assumed that both the observ- 
ables and the wavefunctions are real as happens in most 
cases. As we have seen, for a small quench around criti- 
cality both c„ and p n will be rapidly decreasing after their 
maximal value (in modulus), and a good approximation 
to Eqns. ([5|) and ((6|) can be obtained by retaining only 
few terms. We have observed that the following minimal 
prescription retaining only the three largest components 
works fairly well: 

F(t) = F + Aco8 (w A t)+B cos (uj B t). (7) 

For instance A — 2p pi, B = 2p$p 2 for the Loschmidt 
echo while A — 20o,iCoCi, B — 20a^co c 2 for a more 
generic observable O. The distribution function re- 
lated to the time-signal Eq. 0, P (/) = 8 (f -F{t)), 
has been computed exactly in Ref. [a]. P (/) is a 
symmetric function around the mean F supported in 
\F- P| + \B\\ ,F + \\A\ + \B\\] with logarithmic diver- 
gences at f = F±\\A\ - \B\\ (see Fig.E). 

This scenario can be summarized as follows: For small 
quench around a critical point, generic observables equi- 
librate only very poorly. The distribution function for 
a generic observable is a double peaked distribution with 
a relatively large mean, a behavior completely different 
from the Gaussian one. 

To complete the analysis let us now discuss the case of 
a small quench in a regular point of the phase diagram. 
At regular points there are no gapless excitations and the 
weights are bounded by p (E) < M (E) /A 2 where A is 
the smallest gap. Since the theory is not scale-invariant 
M (E) will not be an homogeneous function, and in par- 
ticular will not display any singularity. The picture then 
is the following: In the perturbative regime (SX 2 L d < 1) 
we still have a "large" lowest weight, but beside po no 
other p n dominates and the sum rule ^ n p n — 1 is satu- 
rated only recurring to a relatively large bunch of p n s. 

In general predicting the precise behavior of observ- 
ables in this case will be difficult as one needs to have 
knowledge of many different weights in Eqns. ([5]) and 
However we can give a simple argument to ex- 
pect a Gaussian behavior for generic case. As we 
have argued, the sum in Eq. ^ contains now many 
terms. If the energy differences E n — E m are rationally 
independent, along the time evolution, each variable 
X n ,m '■= 2(n\0\m)c m Cn will span uniformly the inter- 
val {—2\(n\0\m)c m c n \,2\(n\0\m)c m c n \]. As long as the 
variables X nyTn can be considered independent, O (t) — O 
can be thought of as a sum of independent random vari- 
ables. Since, as we have seen, the sum is made over many 
variables, the central limit theorem applies and the re- 
sulting distribution will be Gaussian. This argument can 
fail when the variables cannot be taken as independent. 
This can happen, for instance, when a certain observable 
is pushed toward its maximum or minimal value by the 
action of some field. Consider for example the case of a 
transverse magnetization crj in presence of a high field 
—hJ2i a i- F° r increasing h the mean of (<r| (t)) will be 
pushed towards one. Since (crj (t)) is supported in [—1, 1] 



P(x| P(m) 




Figure 1: (Color online) Probability distributions for a small 
quench around a critical point. P (x) — S (x — C it)), P (m) = 
5 (m — (af (i))) refer to the Loschmidt echo (upper panel) and 
magnetization respectively (lower panel). The thick lines are 
obtained using the prediction of Eq. (0 using only the three 
largest weights. Note the large spread of the distributions 
compared with their total support: P (x) £ [0, 1] and P (m) £ 
[-1, 1]. Parameters L = 16, ki = k 2 = 0.4, hi = 0.218, Sh = 
h2 — h\ — 0.04, are close to criticality (see HH). The data 
are obtained by Lanczos diagonalization of Eq. |(SJ keeping as 
many lowest energy vectors until the sum rule X^n=o x_1 P" — 
1 was satisfied within prescribed accuracy. 

P(x) P(m) 
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Figure 2: (Color online) Probability distributions for a small 
quench around a regular point. P (x) — S [x — C (t)), P (m) = 
S (m — (a'{ (t))) refer to the Loschmidt echo (upper panel) and 
magnetization respectively (lower panel). The thick lines are 
Gaussian with the same mean and standard deviation. The 
quench is performed in the paramagnetic phase, parameters 
are L — 12, ki = «2 = 0.3, hi = 1.4, Sh = 0.04. Note the very 
small spread of the distributions. The data are obtained by 
full diagonalization of the Hamiltonian Eq. @. Histograms 
are obtained by sampling 40000 random times in an interval 
t £ [0, T] with T = 16000. 



the corresponding distribution can cease to be Gaussian 
as its mean is pushed against the (upper) border of its 
support. In this case the distribution function will look 
like a "squeezed" Gaussian. A similar effect has been 
observed to take place to the Loschmidt echo in Ref. @ 
when the system size becomes the largest scale of the 
system. In any case, however, if the variables cannot 
be considered as independent, any possible distribution 
function (and not only a squeezed Gaussian) can arise. 

Numerical test We will now check our predictions on 
the hand of a non-integrable model. As a test model we 
chose to use the so called TAM Hamiltonian (transverse 
axial next-nearest-neighbor Ising model). The Hamilto- 
nian is 

L 

H = -H («+i - ™* + hof) , (8) 

i=l 



4 



and periodic boundary conditions are used (cr£ +i = of). 
A positive K frustrates the order in the a x direction. The 
reason for our choice is, at least, twofold: i) The TAM 
is a non-integrable generalization of the one-dimensional 
quantum Ising model for which results are already avail- 
able [|[ . ii) The model Eq. ([8|) has only a discrete Z2 sym- 
metry (P z = Y\i of), consequently the ground state lives 
in a large das — 2 L ~ 1 dimensional space. In practice des 
is the effective Hilbert space dimension, and we would 
like it to be as large as possible. For instance, after a 
quench the purity of the equilibrium state is bounded by 
tr(p 2 )> 

^gs- This i s to b e contrasted with other mod- 
els used in the literature with larger symmetry groups 
(i.e. SU (2)) for which the dimension of the block contain- 
ing the ground state is still exponential in L but consid- 
erably reduced with respect to to that of the full Hilbert 
space 2 L . 

The model Eq. (|8|) displays 4 phases (see for instance 
[2lL [22l . [23l . [2H and references therein), ferromagnetic 

+ + ++, antiphase +H , paramagnetic, and a floating 

phase with algebraically decaying spin correlations. In 
particular, for small frustration k < 1/2, increasing the 
external field h there is a transition from ferromagnetic 
to paramagnetic. This transition is believed to fall in 
the Ising universality class, and so the critical theory is 
described by a conformal field theory with central charge 
c = 1/2 and d = Q = v = 1. Wc performed our numerical 
simulation for the critical quench on this critical line. 

We will illustrate our findings for two particular yet 
physically well motivated observables; the Loschmidt 
echo C(t) = J {ijjo\e~ ltH2 \ipo) | and the transverse mag- 
netization m (t) = (ipo (t) I of I ipo (t)). 

Since d — £ = v = 1, according to Eq. ([3]), we ex- 
pect a strong divergence at low energy: p (E) ~ E~ 2 . 
Consequently we expect very few p n , n > to have non- 
negligible weight, and so Eq. ([7]) to be a valid approxima- 
tion. Indeed the results based on numerical diagonaliza- 
tion compare well with the prediction based on Eq. ([7]) 



(Fig. [1} . Note the very large spread of the distributions 
compared to their total support: "poor equilibration". 

For comparison we performed similar numerical simu- 
lation for a small quench in a regular point of the phase 
diagram. As expected the resulting distribution func- 
tions are approximately Gaussian (Fig. [5]). Note the very 
small variances of the distributions already for a rela- 
tively short size: "good equilibration" 

Conclusions In this paper we investigated the de- 
tailed structure of equilibration after a small quench, 
i.e. the system is initialized in the ground state of a 
given Hamiltonian Hi and then let evolve with a slightly 
perturbed Hamiltonian H2 = Hi + 5 XV. In the limit 
5X — > equilibration is trivial in that for all observ- 
ables P (o) = 8{o-{0 (t))f = 5{o-{0)). However this 
limit is approached very differently depending on whether 
Hamiltonian Hi is critical or not. For quenches around 
a regular point of the phase diagram the expected distri- 
bution for generic observables is a Gaussian one. Equili- 
bration arises in the most standard fashion. Instead for 
small quenches around a critical point the situation is 
radically different. The distribution function for generic 
observables P (o) tends to universal double-peaked func- 
tion for relevant perturbations. 

The key step to obtain these results is to character- 
ize the overlaps c„ = (n\ipo) between the initial state 
\4>o) and quenched Hamiltonian eigenstates \n). We have 
shown that, at criticality, the function c(E n ) = c n (E n 
eigenenergy) decays very rapidly: c(E) ~ E^ 1 /^^ and 
this in turns gcncrically implies the observed double- 
peaked distributions. The analytical predictions have 
been checked numerically on the hand of a non-integrable 
extension of the quantum Ising model. 
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